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AN  INVESTIGATION  OF  THE  EFFECT  OF  TWIST  ON  THE  COUPLED 
BENDING-TORSION  VIBRATIONS  OF  A  CANTILEVER  BEAM 


SUMMARY 

An  experimental  study  was  made  of  the  effect  of  twist  on  the  natural 
frequencies,  coupled  in  bending  and  torsion,  of  a  cantilever  beam.  The  shape  of 
the  beam  was  selected  to  provide  a  high  degree  of  coupling.  An  analytic  procedure 
for  determining  these  frequencies  was  also  developed  and  is  included  in  this  report. 

It  was  determined  that  coupling  without  twist  reduced  all  frequencies  from 
their  uncoupled  values.  The  fundamental  frequency  was  found  to  remain  relatively 
constant  as  the  beam  was  twisted.  All  natural  frequencies  higher  than  the 
fundamental  were  lowered  as  the  total  twist  was  increased  to  15  degrees,  but 
remained  relatively  constant  with  further  twist. 


INTRODUCTION 

With  the  advance  of  helicopters  and  turbo-machinery,  the  problem  of  twist 
and  its  effect  upon  the  natural  vibrations  of  a  beam  becomes  of  more  and  more 
interest.  Much  work  has  been  done  to  determine  the  effect  of  twist  on  beards 
having  no  coupling  between  torsion  and  bending,  and  the  effect  of  a  high  degree  of 
bending-torsion  coupling  on  untwisted  beams.   This  report  shows  primarily  the 
effect  of  twist  on  the  natural  frequencies  of  a  particular  beam  highly  coupled  in 
bending  and  torsion.  It  also  shows  the  node  lines  for  all  the  natural  frequencies, 
and  indicates  the  effect  of  bending-torsion  coupling  on  the  natural  frequencies  in 
the  untwisted  beam. 

A  beam  with  dimensions  as  shown  in  Fig.  1  is  considered.  This  particular 
shape  was  chosen  to  give  a  maximum  amount  of  bending-torsion  coupling, 
relatively  easy  tooling  problems,  and  readily  obtainable  sectional  and  elastic 
properties.  In  this  report,  the  term  coupling  refers  to  the  interaction  of  the 
torsional,  flapwise  bending,  and  chordwise  bending  types  of  vibrations.  It  was 
planned  to  use  both  experimental  tests  and  analytic  means  to  determine  the 
natural  frequencies  resulting  from  the  combination  of  twist  and  coupling.  Only 
experimental  results  however,  are  contained  in  this  report.  Difficulties  en- 
countered in  preparing  the  digital  computer  program  prevented  the  inclusion  of 
analytic  results.   The  mathematical  method  used  is  outlined  in  the  Appendix, 

This  method  is  based  on  one  being  used  in  a  report  of  Isakson  and  Eisley 
currently  being  prepared  for  the  University  of  Michigan  Research  Institute.  It 
involves  the  introduction  of  coupling  terms  from  the  equations  of  motion  developed 


by  Houbolt  and  Brooks,  Ref.  1,  into  the  matrix  method  of  solution  of  Targoff,  Ref.2  &3, 
The  Targoff  method  was  described  and  used  in  a  report  of  Isakson  and  Eisley,  Ref.  4. 

The  coordinate  system  used  is  shown  in  Fig,  2.  The  flexural  center  of  the 
cross  section  was  located  according  to  a  formula  from  Roark.  Ref.  5,  which  formula 
was  developed  from  Ref.  7.  The  elastic  axis  of  the  beam  was  maintained  as  a 
straight  line  by  mechanically  twisting  the  beam  about  this  axis.  The  nomenclature 
and  sign  convention  for  the  displacements,  shears  and  bending  moments  used  in  the 
analysis  are  shown  in  Fig.  2.  All  parameters  used  or  needed  are  listed  in  the  table 
of  symbols.  The  non-dimensional  forms  of  these  parameters  used  in  the  analysis 
are  also  listed  in  the  table  of  symbols,  The  mass  and  stiffness  distributions  are 
uniform  for  the  length  of  the  beam. 

In  gathering  the  experimental  data,  tests  were  made  at  total  angles  of  twist 
of  0°,   5.5°,   7.90,   12.9°,  20.7°,  26.7°,   30.4°,  and  40,8°. 

In  conducting  this  investigation,  the  assistance  and  guidance  of  Professors 
Eisley  and  Isakson,  of  the  Department  of  Aeronautical  Engineering  of  the  University 
of  Michigan,  were  of  great  value. 
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EI  ,  EI  bending  stiffness  about  major  and  minor 
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Bending  moment  in  flapwise  direction 
Bending  moment  in  chordwise  direction 
Torque  about  elastic  axis  at  any  cross  section 
Beam  length 
Shear  load  in  flapwise  and  chordwise  direction 


respectively 
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att  mode  constant  lor  determination  of  frequency 

eA  distance  between  area  centroid  of  tensile  member 

and  elastic  axis,  positive  for  centroid  forward 


g  acceleration  due  to  gravity 

k  polar  radius  of  gyration  of  cross-sectional  area, 
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elastic  axis 
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total  loadings  per  unit  length  in     x  J  r\  ,   avid  j 

directions 

total  torque  loadings  per  unit  length  in    X  ,  l\  t  <\n  d     t 

directions 

thickness  of  cross  section  at  any  chordwise  station 

coordinate  in  direction  of    X    axis  measured  from 

the  root  along  the  elastic  axis 

blade  angle  between  major  axis  of  cross-section 

and  the  "y"  axis  for  <fi   -  0    , 

positive  when  leading  edge  is  up. 

rate  of  twist  of  blade  per  unit  length 

increment  in  (j   between  blade  segments 


displacement  of  elastic  axis  in  "y"  aJid  "  z  »  direction, 
respectively 

displacement  of  elastic  axis  in  direction  of  minor  and 
major  principal  axes  of  the  cross -section,  respectively 
total  twist  in  blade  between    X   -  C  and    X  c  K 

cross -sectional  coordinate;  Y\    axis  lies  along  major  axis; 
F   axis  is  perpendicular  to  major  axis  and  passes  through 
elastic  axis. 

E  L, 

angle  of  twisting  deformation,  positive  when  leading  edge  is  up 

mass  per  unit  length  of  blade 

f. 


Note:  all  zero  subscripts  refer  to  root  section 
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Parameter  Values  for  Beam  Under  Consideration 

Values  Independent  of  Twist 

6 
B  93.941  in 

B 

10.5x10"  lb/in2 

EI  in2 

EI  in2 

G                      4.0x10°  lb/in2 


1 
2 

[2e 

n 


% 


93.941 

-11.50 

10.5  x  10 

.7959  x  10 

32.81  x  10 

4.  Ox  106 

1.1864  x  10 

.0758 

3.125 

15.552 

.18723  x  10~4 

1.1864  x  10"  lb   in2 


.  4 
in 


.  4 
in 

.   4 
in 


2  /    2 
lb  sec  /in 

38. 4099  x  10"4  lb  sec2/in2 


J  .2966  in4 


R  26.0  in 


2.350  in 

2 


g  32.2  ft/sec 


k  2.635  in 
A 

1  2.6  in 

m  .  00145  lb  sec2/in 

p  .  00055575  lb  sec  /in2 


Values  which  are  dependent  on  twist 


Total  twist: 

0 

15 

30 

A    (lb  in2) 

1.1864  x 

io6 

1.2699  x 

io6 

1.5196  x  10 

C    (radians) 

0 

-.874780 

-1.396937 

d  (degrees) 

0 

15 

30 

9  (radians) 

0 

.2388 

.4776 

&$  (radians) 

0 

. 02388 

.04776 

/S'(rad/in) 

0 

.00920 

.01838 

Non  Dimensional  Parameters 


Values  independent  of  twist 


498356 


x  10 


-4 


If 


102.2392 


-4 
x  10 


.09038 


.1013 


.1 


? 


For   0 


.0243 


Values  dependent  on  twist 


15 


30 


1.4906 


1.59550 


1.90440 


-.374780 


-1.39694 


12617 


50469 


0472554 


192792 


EQUIPMENT  AND  PROCEDURE 

The  model  used  in  the  present  investigation  was  made  of  2024  T  4  Aluminum, 
machined  from  6"  x  1"  stock.  It  was  machined  to  within  .  015  inch  of  the  given 
beam  dimensions. 

Four  inches  at  the  end  of  the  model  was  left  rectangular  for  mounting  purposes 
(Fig.  1). 

The  bar  was  mounted  with  four  bolts  between  two  modified  channel  sections, 
which  in  turn  were  fastened  to  a  double  box  section  bolted  to  a  structural  member 
in  the  wall  of  the  building.  The  double  box  section  was  made  by  welding  1/2"  inch 
steel  plates  between  the  outer  edges  of  the  flanges  of  a  6  x  6,  5/8  inch  steel  I  beam. 
(Figs.  3  and  4).  The  channel  sections  were  later  reinforced  by  welding  a  0.  5  inch 
steel  gusset  as  shown  in  Fig.  3.  This  was  done  to  observe  the  effect  upon  the 
natural  frequencies  of  substantially  stiffening  the  mounting.  No  measurable  change 
in  frequencies  was  noted.  This  provided  some  indication  that  the  mounting  was 
sufficiently  rigid  to  give  results  that  closely  approached  a  true  cantilever. 

MB  Vibration  Test  Equipment,  Model  Tl  -  32034  was  used  to  vibrate  the 
beam.  This  machine  has  a  frequency  range  of  from  two  to  20,  000  cycles  per  second. 
The  vibrator  itself  was  mounted  to  the  floor  and  fastened  to  the  beam  with  an 
aluminum  tube.  This  tube  was  fastened  to  the  vibrator  with  a  flexure  joint  and  to 
the  beam  with  a  universal  joint  (Figs.  3  and  4). 

Two  universal  joints  were  used  originally,  one  at  each  end  of  the  aluminum 
tube,  but  this  arrangement  allowed  too  much  slack  in  the  connection.  Two  flexure 
joints  were  also  tried  in  the  same  manner  as  the  universals,  but  this  arrangement 


was  too  stiff  to  permit  the  beam  sufficient  freedom  in  torsion. 

To  obtain  the  desired  twist,  special  fittings  were  made  to  fasten  to  each  end 
of  the  bar  so  that  it  could  be  twisted  about  its  elastic  axis.  A  Riehle  torsion  testing 
machine  of  120,  000  in.  lbs.  capacity  was  used  to  twist  the  beam.  The  twisting  rig 
as  shown  in  Figs.  5  and  6  was  used  to  transmit  the  torque. 

The  method  of  determining  natural  frequencies  consisted  of  the  following  steps: 

1)  A  frequency  was  approximately  located  by  listening  for  and  observing  the 
amplitude  maximum.  In  the  case  of  the  fundamental  frequency  this  method  provided 

an  accurate  answer  by  noting  the  frequenoy  at  which  the  measured  maximum  amplitude 
of  the  free  end  of  the  beam  occurred. 

2)  The  beam  had  been  designed  so  that  uncoupled  natural  frequencies  of  torsion 

and  chordwise  bending  occurred  in  the  vicinity  of  a  natural  uncoupled  flapwise 

bending  frequency.  When  the  beam  was  being  vibrated  in  the  vicinity  of  a  band  of 

natural  frequencies,  fine  sand  of  silicone  flour  was  scattered  over  the  surface  of  the 

beam.  The  frequency  was  then  adjusted  until  sharp  node  lines  appeared  and  concided 

with  maximum  amplitudes  of  vibration.  The  maximum  amplitudes  were  recognized 

by  first  turning  down  the  amplitude  control  of  the  vibrator  until  the  sand  or  powder 

was  barely  agitated,  then  adjusting  the  frequency  until  the  maximum  agitation  of  the 

particles  appeared.  This  method  permitted  the  determination  of  amplitude  maxima  & 

in  a  band  of  frequencies  where  sharp  node  lines  did  not  disappear,  but  only  shifted 

their  positions.  Where  several  maxima  ^  appeared  in  a  band,  the  strongest  was 

assumed  to  be  the  frequency  of  the  flapwise  bending  mode.  This  assumption  was 

made  considering  that  a  flapwise  node  was  designed  to  occur  in  eaoh  band,  and  that 

due  to  the  shape  of  the  cross  section,  the  flapwise  vibration  would  be  the  dominant  one. 
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It  is  the  frequency  change  due  to  coupling  effects  as  influenced  by  twist 
rather  than  the  absolute  value  of  the  frequencies  that  is  important  in  this  report. 
Several  factors  occur  in  the  experimental  procedure  which  would  tend  to  produce 
an  experimental  result  different  from  one  arrived  at  by  calculations.  Throughout 
the  experimental  runs,  attention  was  directed  toward  keeping  these  factors  constant. 

The  connecting  rig  between  the  beam  and  the  shaker  added  end  mass  effect; 
this  is  discussed  in  the  next  section  of  this  report.  It  was  noted  that  varying  the 
tightness  of  the  bolts  at  the  mounting  even  slightly  could  produce  a  slight  variation 
in  frequency.  Of  course,  for  all  runs  the  beam  was  mounted  as  securely  as  possible 
so  that  the  net  effect  of  variations  from  one  test  to  another  was  very  small. 

For  the  evaluation  fo  the  analytic  part  of  this  report  the  Royal  Precision 
Electronic  Computer  LGP-30,  manufactured  by  the  Royal  McBee  Corporation  was 
to  be  used.  A  matrix  program  was  set  up  as  shown  in  the  Appendix.  The  computer 
evidently  could  not  handle  the  problem  as  it  was  programmed.  A  great  deal  of 
time  had  been  spent  on  this  part  of  the  investigation  and  insufficient  time  remained 
to  pursue  it  further. 
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RESULTS  AND  DISCUSSION 

The  effects  of  twist  upon  the  natural  frequencies  of  the  beam  are  shown  in 
Figs.  7  through  12.  In  Figs.  7  through  9,  the  absolute  values  of  the  natural 
frequencies  are  plotted  versus  the  angle  of  twist.  In  Figs.  10  through  12,  instead 
of  absolute  frequencies,  the  ratio  of  each  frequency  at  any  degree  of  twist  to  the 
corresponding  value  for  the  untwisted  beam  is  plotted  versus  the  angle  of  twist. 
The  frequency  ratios  are  grouped  according  to  the  assumed  predominant  motion, 
i.e.  flapwise,  chordwise  or  torsional. 

Table  I  lists  the  uncoupled  natural  frequencies  of  the  beam  calculated  on  the 
basis  of  geometry  and  mass  of  the  beam.  These  are  compared  with  the  corresponding 
natural  frequencies  observed  when  vibrating  the  untwisted  beam.  A  coupled  .resonant 
frequency  appears  in  the  frequency  speotrum  corresponding  to  each  calculated  un- 
coupled natural  frequency.  It  was  assumed  that  the  type  of  motion  associated  with 
each  uncoupled  natural  frequency,  remained  as  the  predominant  motion  in  the 
corresponding  coupled  natural  frequency  of  the  untwisted  beam,  and  would  continue 
to  predominate  as  each  frequency  changed  with  twist. 

The  discussion  of  the  effect  of  twist  is  based  primarily  on  Figs.  10  through  12, 
the  frequency  ratio  plots.  In  this  discussion  each  frequency  is  identified  by  its 
predominant  motion. 

Twist  has  very  little  effect  upon  the  fundamental  frequency  (first  flapwise 
bending)  of  the  beam.  This  frequency  increases  very  slightly  as  the  twist  is  increased 
through  30  degrees.  The  data  point  obtained  at  40.4  degrees  of  twist  indicates  that 
there  is  a  possibility  of  a  more  rapid  rise  in  frequency  as  the  twist  is  increased 
beyond  this  point. 
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The  second  flapwise  frequency  drops  steadily  as  twist  is  increased  to 
approximately  15  degrees;  it  remains  constant  as  twist  is  increased  further. 

The  third  and  fourth  flapwise  frequencies  decrease  as  twist  is  increased 
to  10  degrees.  For  the  next  10  degrees  of  twist,  these  frequencies  rise  slightly, 
then  remain  constant  as  twist  is  inoreased  further. 

The  first  torsional  frequency,  corresponding  to  the  fourth  natural  coupled 
frequency  of  the  beam,  drops  off  rather  rapidly  as  twist  is  increased  to  15  degrees, 
it  then  decreases  slowly  as  twist  is  further  increased.  The  second  torsional 
frequency  decreases  until  about  15  degrees  of  twist  is  obtained,  then  tends  to 
remain  constant. 

The  third  torsional  frequency  decreased  to  15  degrees  of  twist  then  appeared 
to  rise  as  twist  was  further  increased.  Data  for  the  third  torsional  frequency  was 
considered  less  reliable  than  the  other  data,  because  of  the  difficulty  encountered 
in  obtaining  its  resonance  point  and  in  obtaining  sharp  node  lines. 

Only  the  first  chordwise  frequency,  corresponding  to  the  third  natural 
coupled  frequency  of  the  beam,  was  obtained.  It  dropped  sharply  until  a  twist  of 
about  20  degrees  was  obtained;  it  then  appeared  to  remain  oonstant  or  rise  slightly 
as  further  twist  was  introduced. 

The  shape  and  location  of  the  node  lines  for  the  different  frequencies  are 

shown  in  Figs.  13  through  20.  These  were  obtained  when  determining  the  natural 

frequencies.  Although  in  some  cases  a  band  of  frequencies  occurred  in  which  the 

node  lines  remained  sharply  defined  and  only  changed  their  orientation  and  position 

with  frequency,  resonant  peaks  within  the  band  were  strong  enough  to  permit 

determination  of  resonant  frequencies.  The  occurence  of  frequency  bands  is 
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indicated  in  Figs.  13  through  20. 

The  shape  and  location  of  the  node  lines  was  significant  in  that  they  helped 
in  the  identification  of  the  type  of  motion.  It  was  noted  that  the  node  lines  tend  to 
move  toward  the  free  end  of  the  beam  as  twist  is  increased  to  about  20  degrees. 
With  further  twist  they  remain  relatively  constant  or  even  move  slightly  away 
from  the  free  end.  The  position  of  the  node  lines  in  each  group  of  frequencies 
reflects  the  flapwise  bending  mode  occurring  in  that  group. 

The  node  line  illustrations  show  the  appearance  of  what  was  assumed  to  be 
resonant  frequencies  of  plate  bending.  This  assumption  was  made  since  no 
natural  frequencies  based  on  flapwise  bending,  chordwise  bending,  or  torsion,  were 
predicted  to  oocur  in  this  range.  This  plate  bending  first  appeared  near  300  cycles 
per  second  and  became  more  pronounced  as  twist  was  increased.  At  high  angles  of 
twist,  (20  to  30  degrees),  it  began  at  approximately  276  cycles  and  produced 
distinct  node  lines  through  a  band  of  400  cycles  per  second.  Twist  had  little  effect 
on  the  mean  frequency  of  this  mode  as  can  be  seen  in  Fig.  21. 

This  plate  bending  appeared  again  at  about  680  cycles  for  a  twist  of  25.  7  degrees, 
becoming  evry  distinct  at  a  twist  of  30.4  degrees.  Another  plate  bending  mode  was 
seen  to  oocur  at  960  cycles  and  first  appeared  at  a  twist  of  5.  5  degrees. 

Fortunately,  the  plate  bending  resonant  frequencies  appeared  between  the 
other  bands  of  frequencies  and  though  they  undoubtedly  influenced  the  shape  of  the 
node  lines  in  the  upper  range  of  frequencies,  it  is  felt  that  they  did  not  greatly 
hinder  the  accurate  determination  of  resonant  frequencies  in  modes  of  vibration 
which  are  of  interest  in  this  paper. 

In  shaking  the  beam,  oertain  frequencies  above  the  fourth  natural  frequency 
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did  not  give  distinct  nodes.  In  these  cases,  the  natural  frequencies  were  plotted 
as  a  band  rather  than  as  a  distintc  point  (Figs.  7  through  9).  This  difficulty  was 
due  not  only  to  the  proximity  of  the  flapwise,  chordwise,  and  torsional  uncoupled 
frequencies,  but  also  to  the  appearance  of  plate  bending  which  was  discussed  above. 
Other  inaccuracies  were  due  to  some  unavoidable  looseness  in  the  vibrating 
mechanism,  primarily  the  universal  joint.  This  looseness  also  limited  the  frequency 
range  obtainable,  since  the  vibrating  movement  was  lost  in  transmittal  at  high 
frequencies.  The  tube  and  fittings  which  connected  the  bean  to  the  vibrator  also 
added  a  slight  end  mass  effect  to  the  beam.  By  adding  additional  known  mass 
weights  to  the  end  of  the  beam,  it  was  determined  that  the  effect  of  the  vibrating  rig 
reduced  the  fundamental  frequency  by  1.  6  cycles  per  second.  The  data  taken  was 
not  corrected  for  this  effect,  however,  since  the  effect  at  higher  frequencies  was 
not  known.  The  frequencies  might  also  have  been  slightly  low  due  to  the  impossibility 
of  obtaining  a  perfectly  rigid  support  at  the  root.  The  base  of  the  beam  however, 
was  considered  sufficiently  rigid  so  that  this  effect  was  minor. 
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CONCLUSIONS 

1.  In  all  cases,  the  effect  of  bending-torsion  coupling,  independent  of  twist, 
appears  to  reduce  the  natural  frequencies  of  the  beam. 

2.  The  fundamental  frequency  is  relatively  insensitive  to  twist,  demonstrating 
only  a  slight  rise  as  the  angle  of  twist  is  increased. 

3.  All  natural  frequencies  above  the  fundamental  one  are  lowered  as  twist  is 
increased  to  approximately  15  degrees. 

4o        Beyond  15  degrees  of  twist  the  natural  frequencies  are  relatively  independent 
of  further  twist. 
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APPENDIX 


METHOD  OF  ANALYSIS 

The  analysis  is  basically  the  one  which  is  outlined  in  Ref.  3  and  the 
Appendix  of  Ref.  4.  It  is  extended  to  include  torsional  vibration  and  bending- 
torsion  coupling  and  is  applied  to  the  non-rotating  case. 

The  beam  is  divided  spanwise  into  ten  equal  segments.  The  mass  of  each 
segment  is  assumed  concentrated  at  its  center,  and  the  bending  stiffnesses  EI 
and  EI    and  angle  of  incidence  B    are  assumed  constant  between  masses.  The 
twist  of  the  beam  is  accounted  for  by  relative  rotations  of  adjacent  uniform  bays 
(between  masses)  about  a  spanwise  axis,  the  change  in  angle,  A  B  ,  being  equal  to 
the  total  twist  in  a  segment  and  occuring  just  outboard  of  the  mass. 

With  torsional  motion  included,  the  following  quantities,  at  any  point  along 
the  beam,  may  be  expressed  as  a  column  matrix: 

Mi 
K 
St 

V2 
M2 

r. 


a     =J 


} 


2 

£2 

Q 

v   0  ; 


The  elements  of  this  matrix  vary  from  one  station  to  the  next  along  the  beam 
in  the  following  fashion:  {A}n  +  |  =  [r]  [e]  [F]  {a}  „  ""ere  [r]  [f£  j  [p] 

are  10  x  10  matrices  representing  linear  relationships,  described  below,  between  the 
elements  of  \&\     and  corresponding  elements  in  <ZMn      •  Coupling  effects  are 
introduced  into  these  relationships  by  consideration  of  the  differential  equations  of 
coupled  motion  developed  by  Houbolt  and  Brooks  in  Ref.  1. 

In  the     £        matrix  which  relates  the  { Af     matrices  across  a  weightless 
section,  the  elements  are  obtained  as  follows: 

(a)   V,<nH>       =       V,(r,) 

wM,w    =   M,(n)  +  \/,Ml 

From  Ref.  1,  at  any  station: 

M,=   EI,  (-v" 5,n/3  +  w" cos fi) 

=   £1,5,"        {f=0      for  eac  h  station  ^ 

(EU")(n^=  n,(n)  +  v,ln)* 
(c,  ^  =         s;M  -  ^l  -  -^v. 

<d)  S,int,)  =    j,ln)-   x,(B»i  +     awf +  y,'"')4 

i     \  ,  x  2. EI,         bETT 

(•>  i4(n+,)  =   vztn) 

n  ntnt,)  =  tip+v^i 

From  Ref.  1,  at  any  station: 

r\isl\  (.""^   -w«*.n/0    -EBip'0' 

=  £  I,  (S':-l(i's;-^  - P'% )  "  E  B, /3'0 


'•*./ 


For  the  untwisted  section     j3  *  0  ,  but  the  second  term  in  the  above  equation 
is  not  dropped  since  it  represents  the  ooupling  between  chordwise  bending 
and  twist  which  can  only  appear  at  this  piont  in  the  relationships.  For  this 
term,  i3   must  equal  the  amount  of  twist  that  actually  exists  in  a  segment 
having  length  x. 

(EI^-EB^'f":    H4W*VaW*  » 

(g)     Q(w,)   =    Qtn) 

From  Ref.  1,  at  any  station: 

Q   =UJ+  EB.H5')1]  $  -   £B4j8'(v*«^  +  w'iin/3) 
=  [  G J  -  E  6,  ( /S'j'J  0'  -  £  B,  J'   ( (5  =0   for  each  stot«n) 

The  #    appears  in  a  ooupling  term  and  is  retained  as  before. 

[GJ  *  £6,(^0'-  £8*  A'  =  QW# 

Equations  (1)  and  (2)  are  solved  simultaneously  for    hi      and      y5       ; 
from     oL       ,      Oi  and       6%       are  obtained  by  integration. 


(2) 


(n) 


f/(n,»)  ^  ^  p'<") 

o2  =         ; -+•         °L 

Elz[/\-  LEM.)Z] 


>(nK  cq  a'r\(nh  *z 


5i      =     ErjA-  SSEH +**'*** 

1  Eli 


0) 


£i* 

3 


where  A  =    GJ    ♦    EB    (/?')* 

Equations  (a)  through  (j)  provide  the  elements  of  the       E       matrix;  these 

elements  are  put  into  non -dimensional  form  and  the  result  is  the  10x10 

[  l]        matrix  whose  terms  are  all  zero  except  for  the  following: 

E-E=E=E=E=E=E=E=E=E  =  1 

11        22        33        44        55        66        77        88        99        1010 

E      =  E       =     I 
21        65 

E=E=E  =        -I 

32        43        87 


E31=  - 

E        =  -1 

41  b 


V 


I 


E42  " 


a2  ll 


E75 
E_      = 


E79 


A-?/ 

76    =  -  1*1     Ijfa) 

ill?    (A    \ 


E85   = 


E„„   » 


86  Z  ft-y 


va-y/ 

^89     **105  2  V    A'f  ) 

106 


E         =  £  tn    /_A_\ 

106         „       ITT/ 

E         =         ^ 

109         4-Y 

The  quantities  B    and  B    are  those  defined  by  Houbolt  and  Brooks,  Ref.  1, 

and  were  evaluated  for  the  particular  beam. 

The  elements  of  the  [fJ  matrix  which  relate  the    <  AJ      matrices  on  either 

side  of  a  concentrated  mass  are  determined  as  follows: 
y/fo")     a      V,(n)    +    m  '£     -I-     m  e*  * 

(b)  m{"h)  =  M,lh)  ♦  avu; 


The  second  term  is  one  resulting  from  the  entire  segment's  undergoing  a 

change  of  slope  when  the  beam  bends. 

k    =     radius  of  gyration  of  a  segment  about  a  chordwise  axis;  m  *  mass  of  a 


segment. 

(o)      S,'""0 

= 

£(»> 

M     V"0 

=8 

6,tn) 

y0>«) 

=7 

v,ln)  + 

W 

ST 

\fcW  4- 

/I/""' 

= 

M1ln)  + 

m  kj  o^ 

(f)  =     Min)  +   rn  kiu)1!^' 

The  second  term  is  a  rotary  inertia  term  as  before; 

k   -  radius  of  gyration  of  a  segment  about  a  vertical  axis;  mcmass  of  a 

segment. 

(g)  dj,  =       bt 

(h)    ^  =   st(n) 

(i)        q("*0    =,      Q<*>  +    me^/n)    +      rn  k/  0H 

where  k    is  the  polar  radius  of  gyration  of 
the  segment. 

0,       0(m,)     =     0(n) 

Equations  (a)  through  (j)  are  put  into  non-dimensional  form  and  provide  the 
elements  of  the  {"fj        matrix: 


F  =F  =F  =F   «F   =  F   *F   =F  »F  =F     =1 
11    22    33    44    55    66    77    88    99    1010 

F.  .*F_.=«  I  P   V 


14    58 

=  F 
110    94 

=  F 
23    67 


F_  ,=F,   =   £  £e„  >* 


The    [h]      matrix  serves  to  rotate  the  coordinate  axes  through  the  angle  A  13 
and  appears  as  follows: 

Rn"R22-R33=R44-R55-R66=R77=R8s'      °°S    ^/3 

R15=R26  =  R37=R48     =     ~  s">  A/3 

R51=R62=R73  =  R84=    ^    A/3 

R99=  R1010    m     X 

The  analysis  then  continues  in  the  same  manner  as  in  Ref .  4,  exoept  that  the 

matrices  are  now  10x10  and  yield  a  10x5  product  when  the  boundary  conditions  at 

the  tip  are  introduced. 

The  final  results  are  obtained  through  evaluation  of  a  5x5  determinate. 

Introduction  of  trial  values  of      Ts    ^        |_    Is    -     r  [u>)    J     ,        produce  points 

which  define  the  residual  curve.  Zero  values  of  this  curve  define  the  natural 

frequencies. 


TABLE  I 
COMPARISON  OF  COUPLED  AND  UNCOUPLED  FREQUENCIES 


Frequency 

Type  of  Vibration 

(Cycles  per  second) 

Uncoupled* 

Coupled 

30.6 

25.9 

1st  Flapwise 

169 

155 

1st  Torsion 

195 

163 

2nd  Flapwise 

198 

177 

1st  Chordwise 

507 

488 

2nd  Torsion 

536 

518 

3rd  Flapwise 

844 



3rd  Torsion 

1054 

925 

4th  Flapwise 

•Does  not  include  the  effect  of  end  masses. 
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